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Background. In 2003, Taiwan experienced a series of outbreaks of severe acute respiratory syndrome (SARS) 
and 1 laboratory-contamination accident. Here we describe a new phylogenetic analytical method to study the 
sources and dissemination paths of SARS-associated coronavirus (SARS-CoV) infections in Taiwan. 

Methods. A phylogenetic analytical tool for combining nucleotide sequences from 6 variable regions of a 
SARS-CoV genome was developed by use of 20 published SARS-CoV sequences; and this method was validated 
by use of 80 published SARS-CoV sequences. Subsequently, this new tool was applied to provide a better under¬ 
standing of the entire complement of Taiwanese SARS-CoV isolates, including 20 previously published and 19 
identified in this study. The epidemiological data were integrated with the results from the phylogenetic tree and 
from the nucleotide-signature pattern. 

Results. The topologies of phylogenetic trees generated by the new and the conventional strategies were similar, 
with the former having better robustness than the latter, especially in comparison with the maximum-likelihood 
trees: the new strategy revealed that during 2003 there were 5 waves of epidemic SARS-CoV infection, which be¬ 
longed to 3 phylogenetic clusters in Taiwan. 

Conclusions. The new strategy is more efficient than its conventional counterparts. The outbreaks of SARS 
in Taiwan originated from multiple sources. 


Severe acute respiratory syndrome (SARS) is caused by 
SARS-associated coronavirus (SARS-CoV) [1-4]. The 
first known outbreak of SARS occurred in China’s 
Guangdong province during November 2002 [5]. By 7 
August of the following year, SARS had spread to >30 
countries, affecting 8096 people and resulting in 774 
deaths worldwide [6]. In Taiwan, the first SARS case 
was diagnosed on 14 March 2003 [7, 8]. This index 
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case involved a Taiwanese businessman who had visited 
Guangdong province during 5-21 February of that year. 
After he returned to Taiwan, he transmitted the disease 
to his wife, his son (SARS-CoV strain TW1), and the 
doctor who treated his son (SARS-CoV strain TW3). 

On 15 March, 7 employees of a Taiwanese construc¬ 
tion company flew from Hong Kong to Beijing; 4 of 
them developed SARS symptoms on 26 March, several 
days after returning to Taiwan [9]. Also on 26 March, 
a man residing at the Amoy Gardens housing complex 
in Hong Kong flew to Taiwan; the following day, he 
took a train from Taipei to Taichung City to visit his 
younger brother. The visitor returned to his Hong Kong 
home on 28 March after having experienced fever dur¬ 
ing the preceding evening. His younger brother (patient 
TWC), who developed symptoms on 31 March, became 
Taiwan’s first SARS-related fatality. 

On 6 April, a Taiwanese woman (patient TW-HP1) 
suffering from fever and coughing that continued for 
several days visited the emergency room (ER) at mu- 
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nicipal hospital HP; she was transferred to another hospital on 
9 April. Seven HP employees, including a laundry worker who 
was identified as the index case, eventually developed SARS, 
resulting, on 24 April, in a shutdown of all operations of hos¬ 
pital HP [7]. In all, 137 probable SARS cases and 26 HP-related 
fatalities resulted from this single nosocomial infection. Patient 
TW-HP1 had not traveled outside Taiwan during the preceding 
12 months, but, on 27 March, she and the visitor from Hong 
Kong had taken the same train; their seats were in separate cars 
(numbers 3 and 5), yet the train ride constitutes their only 
possible point of contact. 

On 28 April, the government of Taiwan imposed mandatory 
quarantines on all air travelers from China, Hong Kong, Sin¬ 
gapore, Macau, and Toronto; however, nosocomial SARS in¬ 
fections continued to be reported in many hospitals islandwide 
[5]. The hospitals that experienced the most-severe outbreaks 
are fisted in figure 1; in the present study, they are referred to 
by initials—“HP,” “JC,” “KC,” “GD,” and “YM.” According to 
Taiwan’s Center for Disease Control (CDC), 346 of the 664 
probable SARS cases that have been reported to the World 
Health Organization (WHO) were confirmed by reverse-tran¬ 
scriptase polymerase chain reaction (RT-PCR) and/or neutral¬ 
izing-antibody tests [10]. Previously, Yeh et al. had studied the 
molecular epidemiology of SARS infection in Taiwan and had 
concluded that the origin of the Taiwanese SARS epidemic was 
mainly either Hong Kong or Guangdong, rather than Beijing 
[11]; in addition, they found that the SARS-CoV isolated from 


the younger brother (i.e., patient TWC) of the visitor from 
Hong Kong was not clustered with other isolates from hospital 
HP [11, 12]. Because a complete genome sequence from the 
nasopharyngeal aspirate from patient TWC (strain TCI) was 
available, we decided to reexamine both (1) the source of in¬ 
fection at hospital HP and (2) the paths of dissemination of 
SARS among hospitals in Taiwan. 

On the evening of 10 December 2003, a medical researcher 
who regularly worked in a biosafety level-4 laboratory started 
to feel feverish [13]. Although he had recently spent 4 days (7- 
10 December) in Singapore, an epidemiological investigation 
indicated that he had contracted SARS from a laboratory-con¬ 
tamination accident on 6 December. According to his descrip¬ 
tion, the SARS-CoV isolates that he had handled in the lab¬ 
oratory included strain HKU-39849 [14] and other clinical 
isolates that have not been well characterized. The nasopha¬ 
ryngeal aspirate from this medical researcher was included in 
the present study. 

The size of the SARS-CoV genome has been measured as 29.7 
kb [15, 16]. A comparative analysis of 14 SARS-CoV isolates has 
identified 2 distinct genotypes, which can be differentiated on 
the basis of 4 single-nucleotide variations (SNVs)—C:C:G:C ver¬ 
sus T:T:T:T—in the variable regions of the SARS-CoV genome 
[17]. Furthermore, the genotype with the T:T:T:T SNVs has been 
associated with infections originating in Hotel M in Hong Kong 
[8, 17]. Currently, phylogenetic analyses of SARS-CoV require 
complete genome sequences [5, 11, 17]; however, because of the 
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Figure 1 . Epidemiological curve of confirmed cases of severe acute respiratory syndrome (SARS) in Taiwan, which Taiwan's Center for Disease 
Control validated by use of either reverse-transcriptase polymerase chain reaction or serological test. Arrows indicate dates of outbreaks of nosocomial 
infection in different hospitals and of diagnoses of SARS in several key patients. 
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limited number of specimens available for complete genome 
sequencing, some researchers have used the SARS-CoV spike 
gene for this purpose; but most results have been less than sat¬ 
isfactory [18, 19]. Therefore, the objectives of the present study 
were (1) to use 20 complete SARS-CoV sequences to combine 
several variable regions of a SARS-CoV for phylogenetic analysis, 
in order to develop a simpler tool; (2) to use a different set of 
sequences from 80 SARS-CoV isolates to validate the method, 
by 3 different phylogenetic analytical methods; and (3) to apply 
this proposed tool to elucidate the origin and paths of dissem¬ 
ination of SARS-CoV infections, on the basis of samples collected 
from hospitals in Taiwan. 

PATIENTS, MATERIALS, AND METHODS 

Patients. Serum, sputum, or throat-swab specimens from 19 
patients with SARS associated with 5 nosocomial infections and 
1 laboratory-contamination incident were collected for the pres¬ 
ent study. For each patient, data on the date at onset of the 
disease and on possible sources of infection were gathered by 
trained interviewers (table 1). We also downloaded, for our 
phylogenetic-tree analyses, 20 complete Taiwanese SARS-CoV 
sequences from GenBank; for 10 of these, SARS-CoV strains 
have been described elsewhere [11]. 

Epidemiological investigation. To track the origins of SARS 
cases in Taiwan, we used SARS-treatment reports written by 
physicians and submitted by their respective hospitals to Tai¬ 
wan’s CDC. We also sent trained interviewers to all hospitals 


that reported nosocomial infections, to gather additional as¬ 
sessment data. We based our investigation on the WHO defi¬ 
nition of SARS [20]. 

Analysis of SARS-CoV sequence variation, and selection of 
variable regions for phylogenetic analysis. Complete nucle¬ 
otide sequences from 20 SARS-CoV isolates available from 
GenBank were aligned by the BioEdit program [21]; the se¬ 
quence variations were analyzed by the SimPlot program (http: 
//sray.med.som.jhmi.edu/RaySoft/SimPlot). Genomic sequences 
from civet-cat SARS-CoV strains SZ-3 and SZ-16 were used as 
standards for comparison [19]. Sequence-variation-distance 
plots were generated by use of an 800-bp window, a 200-bp 
step, and a Jukes-Cantor correction. Initially, 4 variable regions 
(SC18, SC22-23, SC27, and SC28), which contain SNVs of 
different SARS-CoV genotypes [17], were combined for phy¬ 
logenetic analysis. Because the resultant phylogenetic tree was 
not satisfactory, we added another 2 variable regions (SC10 and 
SC20) based on the sequence-variation plot (figure 2). 

Comparison of conventional and proposed strategies for 
phylogenetic analysis. Three methods—neighbor joining (NJ), 
Fitch and Wagner parsimony (Pars), and maximum likelihood 
(ML)—were used for comparison of the conventional (i.e., 
complete genome) and proposed strategies. The proposed strat¬ 
egy for phylogenetic analysis entailed the deletion of conserved 
domains and the combination of minimal variable regions. The 
MEGA2 and Phylip3.6 software packages were used to construct 
the phylogeny [22-24]; and 80 SARS-CoV sequences down- 


Table 1. Demographic data and possible sources of infection in cases of severe acute respiratory syndrome (SARS) that were used 
for molecular epidemiological study. 


Date, 
in 2003, 
of onset 

Patient (sex; age, years) of SARS 


TW-HP1 (F; 47) 
TW-HP2 (M; 41) 
TW-HP3 (F; 37) 
TW-HP4 (M; 31) 
TW-JC2 (F; 38) 
TW-KC1 (M; 54) 
TW-KC3 (F; 42) 
TW-PH1 (M; 60) 
TW-PH2 (F; 51) 
TW-GD1 (F; 26) 
TW-GD2 (M; 75) 
TW-GD3 (M; 73) 
TW-GD4 (M; 75) 
TW-GD5 (M; 78) 
TW-YM1 (F; 47) 
TW-YM2 (F; 67) 
TW-YM3 (F; 90) 
TW-YM4 (M; 86) 
SCVJ (M; 44) 


11 April 

28 April 

29 April 
29 April 

1 May 
15 May 

20 May 
19 May 
23 May 

21 May 
29 May 

2 June 
6 June 
4 June 
8 June 
8 June 
8 June 
8 June 

10 December 


Source of infection 

Transmission during train ride (took same train as visitor from Hong Kong) 

Nosocomial infection (visited emergency room of hospital HP) 

Nosocomial infection (was a patient in hospital HP) 

Family contact (wife was a nurse at hospital HP) 

Nosocomial infection (was radiography technician at hospital JC) 

Family contact (relative was a patient with SARS, at hospital KC) 

Nosocomial infection (was a patient in hospital KC) 

Nosocomial infection (visited hospitals KC and PH) 

Nosocomial infection (provided care to a patient in hospital PH) 

Nosocomial infection (provided care to a patient in hospital GD) 

Nosocomial infection (was a patient in hospital GD) 

Nosocomial infection (was a patient in hospital GD) 

Nosocomial infection (was a patient in hospital GD) 

Nosocomial infection (was a patient in hospital GD) 

Nosocomial infection (provided care to a patient in hospital YM) 

Nosocomial infection (provided care to a patient in hospital YM) 

Nosocomial infection (was a patient taken care of by both TW-YM1 and TW-YM2, in hospital YM) 
Nosocomial infection (was a patient in hospital YM) 

Laboratory-contamination accident 
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Figure 2. Analysis of nucleotide-sequence variation of 20 isolates of human severe acute respiratory syndrome-associated coronavirus (SARS-CoV), by the SimPlot program. Two civet-cat SARS-CoV 
strains were used as references. The X-axis indicates the nucleotide location of the SARS-CoV genome; the Z-axis indicates the rate of nucleotide-sequence variation. Sequence differences are plotted 
with a window of 800 nt and 200-nt steps. The 6 variable regions—SC10, SC18, SC20, SC22-23, SC27, and SC28—used for the combined phylogenetic analyses are in boldface type. 



















































loaded from the GenBank database were used for the com¬ 
parison. A bootstrap analysis of 100 replicates was used to 
compare the robustness of the NJ and Pars trees generated by 
the conventional strategy versus that of the NJ and Pars trees 
generated by the proposed strategy [25]. For the ML method, 
the P value for branch length was calculated by use of the 
hidden Markov model, before the conventional and proposed 
strategies were compared [26]. 

RT-PCR and sequencing. RNA was extracted from serum, 
sputum, or throat-swab specimens by QIAamp viral-RNA mini¬ 
kits (Qiagen). RT-PCR primer pairs for the 6 SARS-CoV var¬ 
iable regions are listed in table 2. RT-PCR was performed in a 
single-tube reaction (Qiagen) using primers from each region. 
The PCR thermocycler program consisted of predenaturing for 
4 min at 95°C; 35 cycles of denaturing for 30 s at 95°C, an¬ 
nealing for 30 s at 52°C, and initial extension at 68°C for 3 
min; and final extension at 68°C for 10 min. The PCR product 
was gel-purified for DNA sequencing by use of an ABI PRISM 
3700 DNA Analyzer (Applied Biosystems). For the laboratory- 
contamination case, multiple primers [17] were used for com¬ 
plete genome sequencing. Superscript III RT (Invitrogen) was 
used for production of cDNA, and a RACE kit (Roche) was 
used to amplify the 5' and 3' ends of cDNA. Both strands of 
the PCR product were analyzed, and the resultant sequences 


were assembled by use of the SeqMan II software package (ver¬ 
sion 5.0; Lasergene). 

Analysis of the nucleotide-signature patterns of different 
waves of epidemic SARS infection. To perform the analysis 
of nucleotide-signature patterns, representative SARS-CoV strains 
with complete sequences were chosen from each wave of ep¬ 
idemic SARS infection in Taiwan. In addition, cases that were 
of unclear origin were also included in the analysis. The nu¬ 
cleotide sequence of SARS-CoV strain Urbani [16] was used 
as a prototype for the comparison. 

RESULTS 

Development of a proposed strategy for phylogenetic analysis 
of SARS-CoV isolates. Pairwise comparisons were used to an¬ 
alyze sequence variations among 20 SARS-CoV isolates. The re¬ 
sults show that the 3' region of the viral genome had the highest 
sequence variation, especially near the junction of replicase lb 
and the spike genes (figure 2). Six variable regions—SC10, SC18, 
SC20, SC22-SC23, SC27, and SC28—were chosen for testing the 
proposed strategy for phylogenetic-tree analysis; the total length 
of the sequences in these 6 regions was 5287 nt. Using 3 phy¬ 
logenetic analytical methods, we compared, for 80 SARS-CoV 
isolates, the topology and robustness of the phylogenetic tree 


Table 2. Primers used for reverse-transcriptase-polymerase chain reaction analysis of 6 variable regions of 
the severe acute respiratory syndrome-associated coronavirus (SARS-CoV) genome. 


Region (np), primer pair 3 

Fragment used 
for combined 
analysis 3 , np (length) 

Coding region 

SC-10 (8964-10100) 

SCI OF: 5'-GGATGCTATGGGCAAACCTGTGCC-3' 

SCI OR: 5'-GGACAGTATACTGTGTCATCCAAC-3' 

9318-9904 (587 bp) 

Replicase 1A 

SC-18 (17002-18124) 

SCI 8F: 5'-CACTCCAAGGACCACCTGGTACTG-3' 

SCI 8R: 5'-CGGTAGGTCATGTCCTTTGGTATG-3' 

17292-17917 (626 bp) 

Replicase 1 B 

SC-20 (18914-20154) 

SC20F: 5'-GGTTGTGAAGTCTGCATTGCTTGC-3' 

SC20R: 5'-CCTCTAAGTCTCTGCTCTGAGTAA-3' 

19111-20047 (937 bp) 

Replicase 1 B 

SC22-23 

SC22 region (20984-22127) 

SC22F: 5-TGACCCTAGGACCAAACATGTGAC-3' 
SC22R: 5'-ACCAGAAGGTAGATCACGAACTAC-3' 
SC23 region (21917-23102) 

SC23F: 5-ACCCATGGGTACACAGACACATAC-3' 
SC23R: 5'-CACACCAGTACCAGTGAGTCCATT-3' 

21237-22830 (1594 bp) 

Replicase 1 B and spike 

SC-27 (26007-27112) 

SC27F: 5'-CAATCGACGGCTCTTCAGGAGTTG-3' 
SC27R: 5-CT 1 CT G CTATT GTAACCTGG AAGT C-3' 

26054-26828 (775 bp) 

ORFs 3 and 4, E and M proteins 

SC-28 (27018-28166) 

SC28F: 5'-AGACCACGCCGGTAGCAACGACAA-3' 
SC28R: 5'-ATGCGGGGGGCACTACGTTGGTTT-3' 

27136-27903 (768 bp) 

ORFs 7-11 


NOTE. np, nucleotide position; ORFs, open reading frames. 
a Nucleotide residues of Tor2 SARS-CoV. 
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Figure 3. Neighbor-joining (NJ) trees of 80 isolates of severe acute respiratory syndrome-associated coronavirus (SARS-CoV) that were downloaded from GenBank, by use of (left) the conventional 
strategy, which uses the complete genome sequence, and (right) the proposed strategy, which uses sequences from 6 SARS-CoV variable regions. Numbers at nodes are bootstrap values (%), and those 
at 3 nodes are circled: node a, between civet-cat SARS-CoV and human SARS-CoV; node b, between SARS-CoV subgroups A1 (early phase) and A2 (middle phase); and node c, between SARS-CoV 
subgroups A and B (late phase) [5], The scale bar indicates the genetic distance, estimated on the basis of Kimura's 2-parameter substitution model [24], 




























































Table 3. Bootstrap and P values for trees generated by use of the complete genome and by use of 6 variable regions 
of sequences of severe acute respiratory syndrome-associated coronavirus, by 3 traditional analytical methods. 


Bootstrap value, % 

Maximum-likelihood 

„ , . Neighbor-joining method Parsimony method method, P 

Sequence(s) used _f_:___ _:_ _ 


for phylogenetic analysis 

Node a 

Node b 

Node c 

Node a 

Node b 

Node c 

Node a 

Node b 

Node c 

Complete genome (29.7 kb) 

100 

100 

63 

100 

100 

56 

<01 

NS 

NS 

6 Variable regions (5.3 kb) 

96 

83 

66 

93 

78 

67 

<01 

<01 

<01 


NOTE. NS, not significant. 


generated by the proposed strategy, which uses only these 6 
regions, versus the topology and robustness of the tree generated 
by the conventional strategy, which uses the complete genome 
sequences. As shown in figure 3, the topology of the NJ tree 
generated by use of the complete genome sequence was almost 
identical to that of the NJ tree generated by use of a combination 
of the 6 variable regions; similar results were obtained for the 
Pars trees and the ML trees (data not shown). 

To compare the robustness of trees generated by the proposed 
method versus that of trees generated by the conventional 
method, we focused on the bootstrap values for 3 bifurcation 
nodes between SARS-CoV clusters: node “a,” between civet- 
cat SARS-CoV and human SARS-CoV; node “b,” between 
SARS-CoV subgroups A1 and A2; and node “c,” between SARS- 
CoV subgroups A and B. In the NJ trees, the bootstrap values 
for nodes a, b, and c were, respectively, 100%, 100%, and 63% 
for the conventional method versus 96%, 83%, and 66% for 
the proposed method (figure 3); similar results were obtained 
for the Pars trees generated by these 2 methods (table 3). For 
the ML tree generated by the proposed method, all P values 
for bifurcation nodes between different clusters were <.01; for 
the ML trees generated by the conventional method, the only 
node with a P value <.01 occurred at the bifurcation node 
between civet-cat SARS-CoV and human SARS-CoV (table 3). 

Phylogenetic analysis of SARS-CoV infections in Taiwan. 
The proposed method is useful for study of the molecular 
epidemiology of SARS-CoV strains, both because of its com¬ 
patibility and because it can obtain results by using fewer spec¬ 
imens; in addition, it is more economical and less time-con¬ 
suming. We applied this tool to better understand the entire 
complement of Taiwanese SARS-CoV isolates, including the 20 
that had been published previously and the additional 19 iden¬ 
tified in the present study. Epidemiological information, na¬ 
sopharyngeal aspirates, and serum samples were collected from 
18 patients with SARS who were treated at 6 Taiwanese hospitals 
(table 1). The SARS-CoV from a laboratory worker (patient 
SCVJ) who contracted SARS during December 2003 was se¬ 
quenced and analyzed. An NJ tree was constructed, and the 
results showed that all 39 Taiwanese SARS-CoV isolates be¬ 


longed to subgroup B and could be divided into 3 clusters— 
Bl, B2, and B3. 

Integration of the results of the phylogenetic analyses and 
the epidemiological information (figure 4) indicated that there 
were 5 waves of epidemic SARS-CoV infection in Taiwan during 
2003. The first wave, during early March, was composed of 1 
imported case, 2 cases of transmission between family members 
(strain TW1), and 1 nosocomial infection (strain TW3); in 
cluster Bl, both SARS-CoV strain TW1 and SARS-CoV strain 
TW3 were clustered with other SARS-CoV strains linked to 
Hotel M in Hong Kong [8, 17]. The second wave (strain TW5) 
consisted of 4 Taiwanese individuals who contracted the disease 
as they flew from Hong Kong to Beijing [9] and who then 
carried it back to Taiwan during mid-March. The third wave, 
which began in late March, consisted of an infection that oc¬ 
curred on a train (patient TW-HP1), a case of transmission 
between family members (strain TCI), and multiple nosoco¬ 
mial infections (all “TW-HP” patients and patients TW-GD1, 
TW-GD-3, and TW-GD4). In cluster B2, all the SARS-CoV 
isolates mentioned above clustered with the SARS-CoV isolates 
from Amoy Gardens (strains CUHK-AG01, CUHK-AG02, and 
CUHK-AG03) [27], with a bootstrap value of 71%. The fourth 
wave, which started during late April and ended in mid-June, 
contained SARS-CoV isolates not only from hospitals JC, KC, 
PH, GD (patients TW-GD2 and TW-GD5), and YM but also 
from sporadic community outbreaks (strains TW10 and TW11), 
and all of these SARS-CoV isolates belonged to cluster B3, with 
a bootstrap value of 71%. The fifth wave occurred during early 
December and began with a laboratory-contamination case, pa¬ 
tient SCVJ. It clustered with both HKU-39849, a SARS-CoV 
strain used in the laboratory [14], and TWC, another strain from 
Taiwan [12]. Sequence-variation rates for the strain from patient 
SCVJ versus strain TWC and for the strain from patient SCVJ 
versus strain HKU-39849 were, respectively, 0.01% (3/29,756) 
and 0.04% (12/29,756). In addition, a 24-nt deletion occurred 
at nucleotide position 26132-26155, resulting in both a frame- 
shift of open reading frame 4 and a deletion of 8 aa residues 
within the small-envelope glycoprotein. 

As shown in figure 5, there were 27 SNVs and 1 dinucleotide 


1484 • JID 2005:191 (1 May) • Lan et al. 


Downloaded from http://jid.oxfordjournals.org/ at Serials Acq (Law) on April 13. 2015 










Phylogenetic clusters 


f 


Subgr. B 


Cluster B3 


Cluster B2 


Cluster B1 


Subgr. A 


TW-YM4 
— TW-PH1 
— TW-KC1 
TW-GD2 
— TW-YM1 

—-*TW-JC2 (4/20/03) 

71 ^vs 

71TW9 
— TW11 
TW10 
TWK 
— TWJ 
TW-KC3 
— TW-GD5 


68 '-. 



. .-HP1 \W\/03) 
_ ' 


TW-HP2 
TW-HP4 

^_m-XS03.r' 

TW3 

B 14 / 031 .^- 

Sin2677 

_ H®39849 




65 r 

Frankfurtl 
Urbani 
TOR2 
Sin2774 

-GZ-C 


Civet Cat [_ 



Infection waves 


Wave 4: Unknown origin. 

Caused nosocomial infections in at 
least 5 hospitals in addition to 
community outbreaks. 


Wave 3: An Amoy Gardens resident 
paid a visit to relatives in Taiwan. 
He was the indirect cause of 
nosocomial infections in 2 
Taiwanese hospitals. 


Wave 2: Four employees of a 
construction company returned to 
Taiwan after a trip to Hong Kong 
and Beijing. 


Wave 1: The index case came back 
from Guangdong province, China. 


Wave 5: A laboratory 
contamination accident. 


Figure 4. Neighbor-joining phylogenetic tree of 39 Taiwanese severe acute respiratory syndrome (SARS)—associated coronavirus isolates, generated 
by the proposed strategy. The 5 waves of epidemic SARS infection that are related to the 3 phylogenetic clusters are described to the right. Numbers 
at nodes are bootstrap values (%). The scale bar indicates genetic distance, estimated on the basis of Kimura's 2-parameter substitution model [24], 


deletion among 14 SARS-CoV isolates; 12 of the 27 SNVs were 
nonsynonymous changes. Distinctive SNVs were identified for 
each wave of epidemic SARS infection in Taiwan—at nucleotide 
position 3165 for wave 1; at nucleotide positions 3852, 11493, 
and 26477 for both wave 3 and wave 4; and at nucleotide 
positions 26203 and 27812 for wave 4. The nucleotide-signature 
pattern between the strain from patient SCVJ (the laboratory- 
contamination case) and strain TWC consisted of 2 SNVs (at 
nucleotide position 16325-26600) and 1 dinucleotide deletion 
(at nucleotide position 27808-27809). 

DISCUSSION 

When, for phylogenetic analysis, we combined sequences from 
different variable regions of the SARS-CoV genome, we assumed 


that no dual infection or recombination between the SARS-CoV 
subgroups had occurred. Because SARS is an acute infectious 
disease, the odds that it will be contracted from 2 subgroups are 
very low. Furthermore, the strategy of combining different ge¬ 
nomic regions for phylogenetic analysis has been used in mo¬ 
lecular epidemiological investigations of the origins of HIV-2. In 
those studies, the evolutionary history of a simian immunode¬ 
ficiency virus/HIV-2 lineage was reconstructed by use of a com¬ 
bination of partial gag and env sequences; the method increased 
the accuracy of the phylogenetic analysis [28]. 

It is noteworthy that the SimPlot analysis demonstrated that 
the 3' region of the viral genome, especially near the junction 
of replicase lb and the spike genes (encoding protein S), had 
the greatest sequence variation (figure 2). The S protein is 
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Figure 5. Signature patterns of single-nucleotide variations in 5 waves of epidemic severe acute respiratory syndrome-associated coronavirus 
(SARS-CoV) infection in Taiwan. Nucleotides are numbered on the basis of the complete genome sequence of SARS-CoV strain Urbani [IB], Deletions 
are denoted by X's. *, Isolate with 24-nt deletion at nucleotide position 26132—26155; **, amino acid residues of SARS-CoV nonstructural protein 
(Nsp) and open reading frame (Orf) [IB]. 


considered to be the most important target for the humoral 
and cellular immune responses to SARS-CoV [15, 16]. 

The 3 algorithms most commonly used in molecular phy¬ 
logenetic analyses are the NJ, ML, and Pars methods; we used 
all 3 to obtain tree topologies. In terms of robustness, a boot¬ 
strap value of 70% is often cited as a cutoff for a reliable cluster 
[25]. Although the bootstrap values for nodes a and b in the 
NJ and Pars trees generated by the conventional method were 
higher than those in the NJ and Pars trees generated by the 
proposed method, the difference was not statistically significant. 
In contrast, the bootstrap values for node c in the Pars trees 
generated by the conventional and proposed methods were 56% 
and 67%, respectively, and the corrresponding bootstrap values 
in the NJ trees were 63% and 66%. This finding was confirmed 
by the ML trees: the P value for node c in the tree generated 
by the conventional method was not statistically significant, 
and that for node c in the tree generated by the proposed 
method was <.01. Accordingly, the tree generated by the pro¬ 
posed method was more reliable; also, because the proposed 
method requires only 7 RT-PCR reactions to perform the analy¬ 
sis, it is less time-consuming and more efficient than the con¬ 
ventional method. To facilitate other laboratories’ future mo¬ 
lecular epidemiological studies of outbreaks of SARS, we have 
made the nucleotide-sequence-alignment file of 80 SARS-CoV 


reference strains available on our center’s Web site (http://www 
.ym.edu.tw/aids/Molepi/). 

In the present study, we used 39 Taiwanese SARS-CoV iso¬ 
lates, including 20 downloaded from the GenBank database, to 
trace the origin and the path of dissemination of SARS-CoV 
infections that occurred in Taiwan during 2003. Phylogenetic 
analyses demonstrated that the Taiwanese SARS-CoV strains 
were distributed in 3 clusters—Bl, B2, and B3—which differ 
from clusters 1-3 reported by Yeh et al. [11]; the latter clusters 
were not defined on the basis of a bootstrap value, whereas we 
used a bootstrap value of 70% as the cutoff to define our clusters 
B1-B3. Therefore, our cluster Bl contains both cluster 1 and 
cluster 2 of Yeh et al., and their cluster 3 was divided into 2 
clusters—B2 and B3—in our study. 

On the basis of the epidemiological data for the patients with 
SARS who were considered in the present study, it is clear that 
Taiwan experienced 5 waves of epidemic SARS infection during 
2003. The first and second waves were in different phylogenetic 
clusters, suggesting that they had different origins. Although 
both the second and the third waves were in the B2 cluster, 
the SARS-CoV isolates in the third wave had their own nucle¬ 
otide-signature pattern (figure 5). Neither the first wave nor 
the second wave led to serious outbreaks, but the third wave— 
originating with a resident of Amoy Gardens who visited Tai- 
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Figure 6. Two possible paths of dissemination—between hospitals on Taiwan Island and a hospital on Penghu Island—of severe acute respiratory 
syndrome-associated coronavirus (SARS-CoV) infection. Numbers in shaded circles indicate the temporal order of nosocomial infections; numbers in 
parentheses indicate dates at onset of the infections. 


wan—led to 1 transmission on a train (strain TWC3), 1 case 
of transmission between family members (strain TCI), and 
nosocomial infections in at least 2 hospitals (HP and GD) in 
northern Taiwan. Only 1 nucleotide difference between strains 
TCI and TWC3 was noted (figure 5). In addition, strain TWC3 
and an Amoy Gardens strain, CUHK-AG01, shared an identical 
sequence, even though the woman who was the source of the 
isolate of strain TWC3 never left Taiwan at any time during 
the epidemic. An epidemiological investigation showed that the 
visitor from Amoy Gardens and this woman sat in different 
cars during the train ride. Because this is the first documented 
case in which there is molecular proof of transmission on a 
train, it raises the question of why only 1 passenger contracted 
the infection. Because it has been reported that SARS appears 
to be most infectious at 6-11 days after onset of illness and 
not during the first day of symptoms [29], we assumed that 
the visitor from Amoy Gardens was not highly infectious during 
the train ride, even though he developed symptoms during the 
same evening that he traveled from Taipei to Taichung. 


It is important to note that TWC, the SARS-CoV strain 
isolated from patient TWC, clustered with WHU, an isolate 
from Wuhan City, China, in cluster 1 but did not cluster with 
either strain TCI (direct sequencing of a sample from patient 
TWC) or strain TWC3 (direct sequencing of a sample from 
patient TW-HP1) (figures 3 and 4). In addition, there was a 
7-nt difference between TWC and TCI (figure 5). If we assume 
that strain CUHK-AG01 represents the first-generation SARS- 
CoV in the transmission link, then both strain TCI and strain 
TWC3 were the second-generation, and TWH was the third- 
generation. According to the SNV-based analysis (figure 5), the 
number of nucleotide changes in the SARS-CoV genome per 
number of intermediate hosts was extremely low (<1 nt 
change/host). Because it has also been shown that no or very 
limited nucleotide changes occur in SARS-CoV sequences from 
either cultures or primary clinical specimens [11, 30], we can 
tentatively conclude that strain TWC is a laboratory contam¬ 
inant and did not originate from patient TCI. 

The origin of the fourth wave (cluster B3) is still unknown. 
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The molecular epidemiological data suggest that it originated 
at hospital JC (patient TW-JC2) and then spread to hospitals 
KC, PH, GD, and YM, as well as to others. The epidemiological 
investigations also support this hypothesis: TW-PH1, a patient 
with SARS who was treated at hospital JC went to Kaohsiung 
City and received treatment at hospital KC, and, after being 
treated at hospital KC (table 1), went back to Penghu Island 
and was hospitalized in hospital PH, where he transmitted the 
disease to other health-care workers. As shown in figure 6, after 
combining the epidemiological data with the results of our phy- 
logenetic-tree analysis, we conclude that the fourth wave of ep¬ 
idemic SARS infection progressed in 2 dissemination paths: the 
first path was from hospital JC to hospital KC in Kaohsiung and 
then to hospital PH on Penghu Island; the second path was from 
hospital JC to hospital GD and then to hospital YM. Two cases 
associated with community outbreaks in Taipei (strains TW10 
and TW11) also belong to this cluster. 

With regard to the laboratory contamination, the laboratory 
researcher said that he had used SARS-CoV strain HKU-39849 
for his experiment; however, our results indicate that the se¬ 
quence in patient SCVJ is more closely related to that of strain 
TWC (figure 5). Because the researcher claimed that he had 
not obtained strain TWC from Taiwan’s CDC and that he had 
used many clinical SARS-CoV strains besides HKU-39849, we 
are continuing our investigation to confirm both (1) whether 
the virus that he had obtained from Taiwan’s CDC was in fact 
strain HKU-39849 or strain TWC and (2) which SARS-CoV 
strain that he handled was the source of the laboratory con¬ 
tamination. Because >30 imported SARS cases have been re¬ 
ported by Taiwan’s CDC but have not yet been analyzed, we 
plan to use the proposed tool to conduct further analyses, in 
an attempt to identify their origins. These contact histories will 
provide valuable information for the control of future SARS 
infections. 

NUCLEOTIDE-SEQUENCE ACCESSION 
NUMBERS 

The SARS-CoV nucleotide sequences (6 sequences for each of 
18 strains) identified during this research have been deposited 
in GenBank (accession numbers AY451856-AY451963). The 
reference sequences (accession numbers) used in our sequence- 
variation analysis were from the following 20 strains: Urbani 
(AY278741), CUHK-W1 (AY278554), TOR2 (AY274119), 
HKU-39849 (AY278491), BJ01 (AY278488), BJ02 (AY278487), 
BJ03 (AY278490), BJ04 (AY279354), GD01 (AY278489), TW1 
(AY291451), TWC (AY321118), SIN2774 (AY283798), SIN2748 
(AY283797), SIN2679 (AY283796), SIN2677 (AY283795), 
SIN2500 (AY283794), HSR1 (AY323977), CUHK-SulO 
(AY282752), Frankfurtl (AY291315), and GZ50 (AY304495). 
In addition, 58 SARS-CoV genomes from GenBank were used 
for comparisons of phylogenetic analyses using the conven¬ 


tional and proposed strategies: TW9 (AY502932), TW8 
(AY502931), TW7 (AY502930), TW6 (AY502929), TW5 
(AY502928), TW4 (AY502927), TW3 (AY502926), TW2 
(AY502925), TW11 (AY502924), TW10 (AY502923), GZ02 
(AY390556), ZS-C (AY395003), LC5 (AY395002), LC4 
(AY395001), LC3 (AY395000), LC2 (AY394999), LC1 
(AY394998), ZS-A (AY394997), ZS-B (AY394996), HSZ-Cc 
(AY394995), HSZ-Bc (AY394994), HZS2-C (AY394992), HZS2- 
Fc (AY394991), HZS2-E (AY394990), HZS2-D (AY394989), 
HZS2-Fb (AY394987), HSZ-Cb (AY394986), HSZ-Bb 
(AY394985), HSZ2-A (AY394983), GZ-C (AY394979), GZ-B 
(AY394978), NS-1 (AY508724), WHU (AY394850), Shang- 
haiQXCl (AY463059), ShanghaiQXC2 (AY463060), GD69 
(AY313906), FRA (AY310120), SoD (AY461660), Sino3-ll 
(AY485278), Sinol-11 (AY485277), CUHK-AG03 (AY345988), 
CUHK-AG02 (AY345987), CUHK-AG01 (AY345986), PUMC03 
(AY357076), PUMC02 (AY357075), PUMC01 (AY350750), 
GZ50 (AY304495), TWC3 (AY362699), TWC2 (AY362698), 
ZMY 1, (AY351680), TWY (AP006561), TWS (AP006560), 
TWK (AP006559), TWJ (AP006558), TWH (AP006557), TC3 
(AY348314), TC2 (AY338175), and TCI (AY338174). Two civet- 
cat SARS-CoV strains (SZ3 [AY304495] and SZ16 [AY304488]) 
were used as the outgroup of the rooted trees [19]. 
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